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^ ' Abstract 

The asymptotic behaviour of a class of inhomogeneous scalar field 
cosmologies with a Liouville type of potential is studied. We define a 
set of new variables for which the phase space of the system of Einstein 
equations is bounded. This allows us to perform a complete analysis 
of the evolution of these cosmologies. We also discuss the extension 
of the cosmic no-hair theorem. 



I Introduction 

Since the proposal by Misner of the "chaotic cosmology program", the 
idea that the Universe emerged from a highly irregular state and that the 
inhomogeneities and the anisotropies were washed away giving place to a 
highly symmetric universe, has been one of the most attractive ideas in cos- 
mology. In spite of the theorem proved by Collins and Hawking which 
states that only a subclass of measure zero of the space of homogeneous so- 
lutions approach isotropy and the alternative approaches that postulated, 
according to the second law of thermodynamics 0, that the universe began 
in a highly regular state, Misner's idea has come back due to the succes of 
the inflationary scenarios. 

The cosmic "no-hair" theorems of Wald |^ for homogeneous models and 
of Jensen and Stein-Schabes for inhomogeneous spacetimes pointed out 
how the introduction of a cosmological constant, which can be considerded as 
induced by a scalar field (inflaton), allows the models to isotropize approach- 
ing the de Sitter solution. This situation, however, may dramatically change 
if one takes into account the dynamical behaviour of the scalar field. Heusler 
0, for example, extending the theorem of Collins and Hawking to the case of 
convex and positive potentials of the scalar field possessing a local minimum 
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has shown that only the homogeneous Bianchi models which admit a FRW 
solution as a particular case approach isotropy. Later Kitada and Maeda ^ 
and Ibanez et al. have shown that when one assumes a Liouville type of 
potential for the scalar field (exponential potential), the Wald's theorem for 
homogeneous solutions still applies if the exponential potential is quite flat. 

Although Jensen and Stein-Schabes extended the Wald theorem to in- 
homogeneous solutions, little is known about the effect of the dynamics of the 
scalar field on the asymptotic behaviour of the models. The first attempts to 
incorporate the effect of the dynamical evolution of the scalar field was made 
by Goldwirth and Piran using numerical studies of inhomogeneous models 
and, later on by Calzetta and Sakellariadou by studying the evolution 
of a family of inhomogeneous Cauchy data. The first inhomogeneous scalar 
field exact solutions with exponential potential of the Einstein field equations 



were obtained by Feinstein and Ibanez [|1T| and it was shown there that the 
scalar field does not guarantee by itself that the model inflates or isotropize. 
Although there are some exact inhomogeneous solutions obtained up to now, 
we are lacking a result similar to Heusler's. In particular, for the exponential 
potential, it would be interesting to study the conditions on the scalar field 
leading to inflation and isotropization. 
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The asymptotic behaviour of homogeneous but anisotropic solutions for 
a perfect fluid has been widely studied. By using the kinematical quantities 
of the fluid as variables, the field equations can be written as an autonomous 



system of ordinary differential equations ||T^. One of the most distinctive 
features is that the equilibrium points of these systems are self-similar solu- 
tions. This analysis has recently been extended to the case when the matter 
source is a scalar field with exponential potential [|13|. Since the scalar field 
is homogeneous, one can globally associate with it a perfect fluid, and then 
by using the kinematical variables of the fluid, the Einstein field equations 
decouple and the phase space of the system becomes bounded. 

The main difficulty in analyzing the asymptotic behaviour of inhomoge- 
neous solutions is related with the fact that the field equations involve partial 
derivatives. The study of such systems was initiated in [Q, where a partic- 
ular family of self-silmilar solutions with perfect fluid, for which the Einstein 
equations reduced to an autonomous system, was studied. Since the source of 
the metric was a perfect fluid, the kinematical quantities could be used again 
to analyse the asymptotic behaviour of the system, as in the homogeneous 



case. 



In dealing with an inhomogeneous scalar field one can not apply the 
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analysis developped in [T^ due to the fact that one can not globally interpret 
the scalar field as a perfect fiuid. Therefore in order to investigate whether 
the inhomogeneous scalar field cosmologies undergo an infiationary epoch 
leading to the homogenization and the isotropization of the spacetime, one 
does have to look for a different way to tackle this problem. 

In this paper we initiate the study of the asymptotic behaviour of the 
scalar field cosmologies by considering, as a first step, a 2 parameter family 
of G2 self-similar solutions. The fact that the scalar field is not equivalent 
to a perfect fiuid prevents the use of the kinematical quantities to describe 
the evolution of the solutions. Despite of the lack of a preferred timelike 
congruence in the spacetime we manage to find a set of new variables, in terms 
of which the phase space of the system becomes bounded. This allows us to 
perform a complete analysis of the asymptotic behaviour of these spacetimes. 
In addition, the way to introduce the set of new variables is a promising 
method for dealing with more general solutions. 

The plan of the paper is as follows: in Section 2 we present the metric 
and introduce the new variables. Section 3 is devoted to the analysis of the 
phase space and in Section 4 we discuss the results. 
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II The metric and the compactified phase 
space 

We will consider solutions with one- dimensional inhomogeneity. These met- 
rics are described by the generalized Einstein-Rosen spacetimes which admit 
an Abelian group of isometries G2- If the two Killing vectors are hypersurface 
orthogonal the line element can be written as 

ds" = (-de + dz^) + G (e'^ dx^ + e"''dy^) , (1) 

where the metric functions depend on t and z and the Killing vectors are 
and dy. 

The matter source for the metric is that of a minimally coupled scalar 
field with potential, for which the stress-energy tensor is given by 

Tab = 0,a0,6 - gab i^^ A'" + ^(0)) , (2) 

(latin indices run from to 3) with the Liouville type of the potential 

y(0)=Ae'='^, A>0. (3) 



It is well known that as long as the gradient of the scalar field remains 



timelike, (2) can be rewritten as a perfect fluid stress-energy tensor. 



To simplify the equations we will concentrate in this paper on the class 
of solutions for which the element of transitivity surface G is homogeneous 

G = G{t), (4) 

which is suitable to a description of cosmological models. 



It was shown in |]T6[ that if one assumes separability of the metric and 
separability, in the additive sense, of the scalar field one obtains, from the 
field equations, that the dependence of these functions on the variable z is 
linear 

h{t,z) =p{t)+az, F{t,z) = f{t) + cz, (j){t, z) = (f){t) + bz, (5) 

where a, b and c are arbitrary constants which drive the inhomogeneity. 
This linear behaviour for the inhomogeneity has been recently considered, in 
a different context, by Vilenkin [T^. The general solution for this class of 



metrics and the study of few particular examples were given in |Tl[ . 

Before going to the field equations it is worth mentioning that the metric 
(l)-(5) admits a homothetic vector field given by 

H=- — + {l--)x — +{l + -)y—. (6) 



c dz \ cj (9x \ cj dy 



It was conjectured in [|14[ that the self-similar solutions could be the at- 
tractors of the G2 cosmologies and in this sense the metric (l)-(5) would play 
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an important role in the study of the asymptotic behaviour of the inhomoge- 
neous solutions. The metric (l)-(5) is similar to the metric studied by Hewitt 
et al (it is obtained interchanging t by 2; in (5)). Both metrics, however, 
differ in the character of the orbits of the similarity group and in the type of 
the source. In our case the orbits of (6) are spacelike given by t =constant 
and the scalar field is inhomogeneous and, therefore, it is not equivalent, in 
general, to a perfect fluid. 

The Einstein equations and the Klein- Gordon equation for the scalar field 
with the ansatz (5) are given by 

^-2V{<P)e^ = 0, (7) 
G 

P+-P = 0, (8) 
G 

ap-c- + 2b^ = 0, (9) 
G 

i + '^i + kV{4,)e'' = 0. (12) 

For a massless scalar field, ^ = 0, the former system of equations turns 
out to be easily integrable giving the following solution 

p{t) = A\nt, (j){t) = B\nt, G{t) = t, 



m = (]A^ + S^-^)lnt+l(la^ + 6^)^^ (13) 



where A and B are constants subject to the condition 

aA-c + 2bB^Q. (14) 

When V 0, from (7) or (12), and due to the exponential form of the 
potential 

c^-kb. (15) 

In this case if the constant 6 = then, from (9), either p — and the metric 
is Bianchi type VIq, or a = and the metric is Bianchi type I. Let's note 
that whatever value V takes, when a = b = the metric is Bianchi type 1. 

In order to look for a suitable set of new variables to compactify the phase 
space we write the generalized Friedman equation which is obtained from (7) 
and (10) 

" K^) ' " ^ ^"^^ ^^'^ ^ ^ ^^^'^^"^ ^ °' ^^^^ 

This equation suggests the introduction of the following set of variables 

P 



P = ^ — J > * = a/2- 



G/G + f' G/G + f' 

* = r = i^. (17) 

G/G + f G/G + f 



In terms of these new variables (16) is written as 

1 _ ^2 _ ^2 _ ^2 _ ^2 ^ + > . (18) 

{gig + /) 

Therefore the space of solutions described by the variables (17) is bounded 
to the inside of the 4-sphere (18). Let's note that, from (16), GjG + / has 
to be different from zero, unless A be negative. 

In dealing with perfect fluid models new variables were defined as the 
kinematical quantities of the fluid divided by an appropiate power of the 
rate of expansion 9. This assures a good behaviour of these variables near 
the initial singularity. Let's note that our new variables (17) are divided 
by the quantity G/G + f which is not related with the expansion of any 
timelike congruence. Nevertheless, from the general behaviour near the ini- 



tial singularity found by Belinskii et al [18| and from the work of Isenberg 



and Moncrief [T^ we can assume that the metric (l)-(5), near the initial 
singularity behaves, for each value of the coordinate like a Kasner model. 
Therefore, when t — * 

G ^ t, /~p~(/)~lnt, (19) 

and the variables (17) remain bounded when t tends to zero. 

From (18) we see that the points on the surface of the 4-sphere represent 
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either homogeneous Bianchi type I solutions (when constants a and b are 
zero) or the initial singularity of the models. 
By using (17), (9) is written as 

0/3 + 6^2^ -c(l -$) = 0. (20) 

Except for the trivial case when a — b — 0, this equation gives a constraint 
for the constants a and b. Alternatively, if a and b are fixed one can look 
at this equation as giving a plane, intersection of which with the sphere (18) 
describes the phase space. In the study of the equilibrium points in the next 
Section a and b will be arbitrary chosen constrained by (20). 
Using the variables (17), (7) becomes 



G 1 / r \^ /gV 



(21) 



G 2 vi-$y \GJ 

This equation decouples from the rest of field equations if we introduce a 
new time coordinate 

^^=^ + / = ^^. (22) 

dt G ■' 1-^ G ^ ^ 

Near the initial singularity 

dr 1 , , 

^ .^-oo. (23) 

thus T varies from — oo to +oo. In terms of this new time the field equations 
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are written as 

(24) 
(25) 
(26) 
(27) 

where ' means derivative with respect the new time r. Differentiating (20) 
with respect to r one can easily see that the constraint equation (20) holds 
for all values of r, as long as the initial conditions verify the equation (20). 
Threrefore (24)- (27) along with the constraint equation (20) describe the 
evolution of the metric (l)-(5). 





-P (i - 


_ ^2 _ ^2^ 






(l - 


_ ^2 _ ^2^ 


' 2x/2 ' 




(1-$) (l- 


_ ^^2 _ ^2 _ 




r' = 


-r - 


^2 _ ^2 _|_ 


1^ ^ T 
_$ ^ 

2 2^2 



III The Equilibrium Points and The Invari- 
ant Sets 

In this Section we shall study the qualitative behaviour of the trajectories 
of the system (24)-(27). Let's first note that the system admits the discrete 
symmetry T —>■ —T, and, therefore, without loss of generality the study of 
the equilibrium points will be restricted to F > 0. 
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Equilibrium points 

The equilibrium points of an autonomous system play an important role in 
describing the qualitative behaviour of its solutions. The local stability of the 
equilibrium points are given by the eigenvalues of the linearized differential 
equations. The equilibrium points of the system (24)- (27) can be found 
explicitly and we now give them and their character. 



2 + A;2 2 + F 2 + Fj 

This point is outside the sphere when fc^ < 2 and it is on the surface of 
the sphere when A;^ = 2. The constraint equation (20) is trivially satisfied for 
all values of the constants a and h. The corresponding metric and the scalar 
field when k"^ > 2 are 

k 

= --At + bz, (28) 

where 



^ = T' "^ = 1/ P-2 • 



This metric was obtained in [|lT| and is of Bianchi type VI. The eigenvalues 
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of the linearized system are 

-2 -2 -1 + V5 - 2fc2 -1 - V5 - 2A;2 



(30) 



F + 2' F + 2' F + 2 ' + 2 

When k"^ < 5/2 the equihbrium point is a stable node and when A;^ > 5/2 a 
stable focus. 



1 2v^ 2 2V2 



This equilibrium point is on the surface of the sphere. It disappears 
when k^ > 6 and coincides with the former equilibrium point when k^ — 2. 
Substituting the former values of the variables of the equilibrium point into 
(17) one easily gets that the solution corresponds to a homogeneous Bianchi 
type I solution (a — b — 0). The line element and the scalar field are 

= -^1-^ (31) 

where 

C ^ 2 (6 - k^) 

A(F-2)" ^ ' 

and the metric represents the FRW universe with massless minimally coupled 
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scalar field as a source. The eigenvalues of the hnearized system are 

-^(6-^^), -^(2-^^), 4^^-^')' ^^^^ 

When < 2 this equilibrium point is a stable node but when 2 < A;^ < 6 is 
a saddle point. 

It is important to note that the constraint equation (20) is satisfied not 
only because a and b vanish but because (3 = and \/2^ + — $) = 0. 
If we look into the evolution in time of a particular solution with a and b 
different from zero, and with A;^ < 2, /3, ^ and $ will take values such that 
the constraint equation will be always satisfied and as t — > oo the solution 
will approache this equilibrium point becoming, therefore, homogeneous and 
isotropic. 

Ill + ^2 ^ $2 ^ 1, r = o} 

This ring of equilibrium points belongs to the surface of the sphere and 
represents, therefore, Bianchi type I solutions with a minimally coupled 
massless scalar field. By choosing a particular point on the surface, i.e. 
/3 = Oi, ^ = 02 and $ = 03 with + a| + a| = 1, the eigenvalues of the 
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linearized system are 



2-03 + ^02 , 2(1-03), 0, 0. 



(34) 



Hence, any of these points are unstable, except when $ = 1 (03 = 1), which 
corresponds to the Minkowski spacetime. 

Invariant sets 

Besides the equilibrium points, the existence of invariant sets help to describe 
the qualitative behaviour of the solutions of an autonomous system. In our 
case, there are three invariant sets, two of them describing massless and 
exponential potential scalar field solutions respectively while the third gives 
the Bianchi type I solutions. 

The points with F = (massless scalar field spacetimes) compose an 
invariant set of the dynamical system (24)-(27) whose solutions are given by 
(13) and (14). The dynamics of the points of this subspace can be easily 
studied: when F = the dynamical system reduces to 





(35) 



(1 - *) (1 - - ^2 _ ^ 
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with straight hnes as solutions 



^0-1 



■p + h 




(36) 



These hnes start on the surface of the sphere + + $^ = 1 and intersect 
at the point p — ^ — 0,^ — 1. Thus the solutions evolve from the Kasner 
initial singularity to the Minkowski spacetime. 

The second invariant set is given by the points with F > (and by 
symmetry F < 0) , with equilibrium points being the points I and II described 
above. 

Finnally the third invariant set is given by the points on the surface of 
the sphere (/3^ + + + F^ = 1). In this case the dynamical system is 
reduced to 



and the equilibrium points of the system (37) are the points II and III. 





(37) 



The solutions of this system are again straight hnes. 




(38) 
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The behaviour of the system (24)-(27) can be visuahzed in Fig.l where 
the phase space for /? = and F > is depicted for different values of 
the constant k. The positions of the equihbrium point I as function of k 
are represented by the dashed hne. All the solutions start on the circle 

+ = 1, r = 0. When F = (massless scalar field), for all the cases, 
the solutions tend to the point $ = 1, = 0. When k"^ < 2 the point I is 
outside the sphere and the solutions tend to the equilibrium point II. When 
2 < A;^ < 6 point I is inside the sphere and the solutions evolve either to the 
point I (those that go inside the sphere) or to the point II (those that he 
on the boundary of the sphere). When k^ > 6 the only attracting point is 
point I (inside the sphere). The trajectories on the surface start and finish 
on the circle <l>^ + = 1, F = in such a way that their projections on the 
plane F = are straight line directed to the point ^ — — = | which is 
outside the surface. 

The behaviour described above remains the same when (3 is different from 
zero. 
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IV Conclusions 

We have studied in this paper the asymptotic behaviour of a particular class 
of inhomogeneous solutions with a minimally coupled scalar field with an 
exponential potential. The metric belongs to the class of G2 cosmologies. 

We have succeeded to define a set of new variables for which the entire 
phase space is bounded by a 4-sphere. From the analysis of the dynamical 
system obtained for this metric, we deduce the following: 



i) As in the homogeneous case |T^, the dynamical behaviour of the met- 



ric depends on the parameter k which is related to the mass of the scalar 
field. The parameters a and h which drive the inhomogeneity do not play a 
significant role in this behaviour. 

ii) When k"^ < 2 the only equilibrium point is that given by the FRW uni- 
verse. The trajectories evolve from the surface of the 4-sphere, representing 
the Kasner regime near the initial singularity, towards the isotropic equilib- 
rium point II. When k"^ > 2 there are two equilibrium points, one is the 
homogeneous Bianchi type VI (point I) and the second is the FRW solution 
which is a saddle point and, therefore, unstable against small changes of the 
initial conditions. This means that the cosmic no-hair theorem holds in this 
case provided k"^ <2. 
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The behaviour described above is almost identical to that of homogeneous 



models |]T3|, indicating that the introduction of the inhomogeneity through 
the constants a and h does not affect remarkably the dynamical system. 
It is likely that the presence of a more "strong" inhomogeneity, assuming 
another dependence on the spatial coordinate, may change the behaviour of 
the solutions and the conditions to isotropize. We consider, therefore, that it 
is important to extend this analysis to a more general class of inhomogeneous 
scalar field solutions. 

As to the inflation of the models, for inhomogeneous solutions there 
is no natural way to see whether the models inflate. The hypersurfaces 
=constant do not define a globally timelike observer and, therefore, one 
should use a weaker way to specify the inflationary behaviour. It has been 



suggested [£0|, for example, to look at the fulfillment of the strong energy 
condition, the breaking of which is a necessary condition for a model to 
inflate. The energy density and the pressure of the scalar field are given by 

P=-\^,ar + V, P=-\^,ar-V. (39) 

In terms of the new variables, the breaking of the energy condition is written 
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as 




G/G + f 



262 




For the equilibrium points I and II the expression 



2{e - 2) 

(2 + A;2)2 



and 



3(P_2) 
16 



(41) 



respectively. That means that if < 2 both points describe spacetimes 
inflating, but if A;^ > 2 the energy condition could be broken only by the so- 
lution represented by the point I depending on the value of b. This behaviour 
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is similar again with that of the homogeneous models Q. 
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Figure Captions 



Figure 1.- Phase space when f3 — and F > 0, and for different values of 
the constant k. The equihbrium point I is represented by a square and the 
equihbrium point II by a circle. The dashed line describes the position of the 
point I for different values of k. Faded lines represent the trajectories of the 
solutions lying on the surface of the sphere which is an invariant set. The 
plane F = is an invariant set describing the massless scalar field solutions. 
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— Figure 1.b — 
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